Skyrmion vibrational energies with a generalized mass term 
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We study various properties of a one parameter mass term for the Skyrme model, originating from 
the works of Kopeliovich, Piette and Zakrzewski through the use of axially symmetric solutions 
obtained numerically by simulated-annealing. These solutions allow us to observe asymptotic be- 
haviors of the B = 2 binding energies that differ to those previously obtained [9| . We also decipher 
the characteristics of three distinct vibrational modes that appear as eigenstates of the vibrational 
Hamiltonian. This analysis further examine the assertion that the one parameter mass term offers 
a better account of baryonic matter than the traditional mass term. 

PACS numbers: 12.39.Dc, ll.lO.Lm 



INTRODUCTION 



CN . Being one of the primary candidates for an effective low-energy theory of QCD, the Skyrme model [l|, 0, 0] has 
relished from a large amount of study after which it was realized that it possessed the same symmetry properties as 
' Cj 'i QCD in the limit of large Nc [5]. The topological solitons that appear as solutions to the model's equations of motion 
O j' are identified as mathematical representations of nuclear matter. These solitons, or Skyrmions, are then quantized to 
obtain physical properties of nuclei. Studies have shown that for the nucleon 5 = 1, these calculated properties are 
within a 30% margin of error from experimental data [6]. 
fH I However, when one raises the baryonic number to any 5 > 2, the solutions that have been obtained thus far do not 
correctly describe the presumed geometric properties of nuclear matter. Experimental data indicates that the nucleons 
seem to preserve their individuality within nuclei. For example, the B = 2 solutions have the geometrical shape of a 
toroid, whereas the deuteron, the only stable nuclei with 5 = 2, has the presumed shape of two deformed nucleons 

■ lightly bound together. Another problem is that the binding energy of the toroidal Skyrmion is much too high, roughly 
^^r) \ ^ 80 MeV (Recall that the deuteron has a binding energy of 2.224 MeV), which is also a likely contributing factor to 
CO ■ the odd toroidal shape it possesses. Even worst, both the geometrical shape and binding energy problems persist and 
CO ' are amplified as we increase the value of B. For these reasons, the solitonic field configurations within the standard 

. Skyrme model can be viewed as too malleable in the sense that Skyrmions deform noticeably as they form a bound 

■ state. For the purpose of comparison here we shall use vibrational energies as a quantifiable measure of rigidity (as 
0^ ! opposed to malleability) of the field and analyze the ratio -^j^^ where i = 6r, 2, 3 labels the vibrational energy of the 

■ breathing mode and the two other eigenmodes, and where E^^'^ is the total energy of the soliton (more on this in Sec. 
III). Reaching for higher rigidity solutions can also be understood intuitively by the fact that every nucleon within a 
nuclei of 5 = 2 or higher should presumably be deformed while maintaining its individuality within the nucleus, and 
not meld with others to form complex geometrical objects. A pragmatic goal towards improving the model would be 
to find extensions to the current model that point toward a solution to these issues. Therefore, in this optic, one may 
consider various generalizations of the original Skyrme Lagrangian to find which types lowers the binding energies. In 
this paper, we limit our analysis to a one-parameter generalization to the standard mass term of the Skyrme model. 

In the following, we begin by reviewing the standard Skyrme model as well as reveal the studied generalized 
mass term that introduces a dimensionless parameter labeled D. Initially considered and analyzed by Piette and 
Zakrewski this new mass term is studied further in the context of rigidity here using the simulated-annealing 
numerical algorithm to accurately and effectively minimize the pion field configuration. These exact solutions allow 
us to deepen our understanding of the dependency between binding energies and the parameter and will also 
enable us to conclude that solitonic solutions obtained through the rational map ansatz Q are not good indicators 
of this dependence (Sec. IV). We consider solutions for 5 = 1 and B = 2 with axially symmetric configurations. For 
5 = 1, axial symmetry is an exact symmetry both for the static solution and the rotationally deformed solution. Both 
calculation will be performed providing a quantitative insight on such deformations otherwise shown to be significant 
[lO, [O, [HI- For 5 = 2, the situation is somewhat different. Rotational deformation breaks axial symmetry and 
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would in principle require a full 3D computation. So only the static energy will be minimized for the B = 2 case even 
though axially symmetric solutions were found to be a good approximation in [12]. This was done in order to make 
the numerical efforts more tractable and time efficient. Furthermore, the same solutions are used to compute how the 
vibrational energies of three eigenstates behave as we increase D. The methods used for this analysis is outlined in 
Sec. Ill, and the results indicate quantitatively that the mass term increases the rigidity of Skyrmions with B = 1 
while its effect on the binding energy depend on (iso-)spin. 



II. THE SKYRME MODEL WITH AN EXTENDED MASS TERM 

The model initially proposed by Skyrme comprised of the two term Lagrangian 



L= I (fx 



(1) 



where = U^d^U are the chiral currents associated to the three component pion fields tt such that the SU (2) matrix 
U is defined by 

U = a ^ir ■ TT. (2) 

and e are respectively the pion decay constant and the dimensionless Skyrme parameter. The r's found in ([2]) are 
simply the Pauli matrices, and the field a is an additional scalar field that must satisfy the constraint + tt • tt = 1 
in order to avoid adding unphysical degrees of freedom and to enable the possibility of having solitonic solutions. 
Each solution to ([1]) with boundary condition 

U{r,t)^l as |r| ^ oo (3) 

fall into distinct topological sectors that are distinguished by their given topological charge, or baryon number. For 
5 = 1, spherical symmetry can be utilized to extract a convincing picture of a nucleon through the use of the ansatz 

JJ _ ^iTiXiF{r) ^ 

where one only has to solve an ordinary differential equation involving the chiral angle F{r) that ([T]) brings about 0. 

However, the physics that ([Tj) describes lacks in many respects. First, it is possible to add as many higher order 
terms of the form Tr ([L, L]^) as one wishes, to account for all possible interactions, and it is simply unknown how the 
sum of such terms would affect the resulting solutions although it is generally assumed that higher order terms 
in derivatives could be neglected in the low-energy limit. One may also wonder if the two terms of ([T]), among all the 
possible terms, are really the two optimal ones for portraying the physics of nuclear matter. Second, all solutions for 
B > 2 obtained thus far do not correctly characterize the presumed geometric properties of nuclei that are found in 
nature. Third, it does not take into account the mass the of pions, but fortunately, this problem is easier to resolve. 
One can simply add a chiral symmetry breaking term proportional to the pion mass squared, which has first been 
successfully introduced by Adkins and Nappi [7], in the form 

^^^Tr (1-/7), (4) 



where m-j^ is the pion mass (We assume that the three pions tt^, 7r+, and 7t~ are of equal masses). This term has had 
the added benefit of eliminating shell-like configurations and favoring energy densities that are higher at their centers, 
which is more appealing since it is known that nucleons have roughly even matter densities within their shell radii. 
Hence, the Lagrangian ([1]) together with the mass term @ is what we consider to be the standard Skyrme model. 
Of course, as we have mentioned, many additions and modifications can be made, particularly to the mass term (|4|). 

If we only constrain ourselves to mass terms that obey the boundary condition (|3|), Kopeliovich, Piette and Za- 
krzewski [8] have shown that a generalized mass term can have the form 

-^Tr^l-l^MU'*). (5) 

where the function g{p) and constant K must obey 

/ 9{p)dp = 1 and K = g{jp)p^dp. (6) 

^ oo J oo 
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However, of the many mass terms that are evidently possible, we will study a particular one-parameter family that 
has the additional property of disfavoring shell-like configurations as did the standard mass term (jH). One then hope 
that this new term might also improve the overall properties of Skyrmions, such as providing a better account of their 
binding energies. This one-parameter mass term is 



(1 - 5D) 



Tr(l -U-D{U^ - U^)), 



(7) 



is based on the function 



l{p) = 5{p-l) + D {S{p - 2) - 5{p - 3)) 



where K = 1 — 5D and the parameter D can span the range [0,0.2 [ . At D = 0.2, we note that the mass term is 
infinite. 

Following their original proposal of the general term ([5]), Piette and Zakrzewski [9] studied ([7j) with the use of 
rational map (RM) solutions. In Sec. IV, we will show the results of a similar study that use exact numerical 
solutions instead of approximate RM map solutions. This will demonstrate the limitations of the rational map ansatz 
when it comes to give a precise measurement of the binding energies as a function of especially as it approaches 
the critical value of 0.2. 

Before we discuss binding energies, we will describe in the following section how the study of vibrational modes, a 
worthy study in its own right, can lead us to quantitatively understand how the B = 1 solutions gradually get more 
rigid as D increases. The methods used in this paper for calculating the vibrational energies were previously outlined 
by Had.juk andt Schwesinger [15], which can be contrasted with the techniques of Barnes et al. [13 and of Lin and 
Piette [iQ] where the latter performed their vibrational calculations using rational map solutions. 

III. VIBRATIONAL MODES AND THEIR ENERGIES 

Originally put forward by Hajduk and Schwesinger [15], the method for obtaining the vibrational modes and their 
energies begins by performing the global scale transformation 



a{x) cr{(3kXk) and 7ra{x) i^aiPk^k), 



(8) 



where the jSk are scaling parameters, and where we have assumed that the scaling is uniform with respect to the 
Cartesian coordinates x, y and z. The advantage of such a scaling lies in its simplicity both in its mathematical 
and numerical treatment. Its disadvantage is that it is certainly not as general as would be an arbitrary local scale 
transformation. However, it allows to compute vibrational modes that probe the global rigidity of Skyrmions, giving 
us a clearer indication whether or not (and how) the field configurations are getting more rigid with respect to the 
parameter D. Substituting ([8|) into ^ together with the new mass term (O yields a Lagrangian of the form 



L 



(9) 



where the matrices Mij{P) and V{P) are obtained by direct inspection after all substitutions are made. Since we are 
only concerned with small amplitude oscillations, we perform an expansion around the minimum of V by taking 



which results in the expansion of Mij and V as 



(10) 

(11) 



Mij{/3)=M, 



{■qkdkMi. 



V(f3) = V 



1 



PI 



PI 

ViVj{didjV) 



/3?,/30,/33° 



+ ... 



r 







1 



(12) 



(13) 
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L = :^M^j{P)mrij - -mr]jV,j. (14) 



Keeping only terms up to order 0{r]^) gives us the Lagrangian 

In order to compute the vibrational Hamiltonian we must now find a coordinate transformation that satisfies 

A^M^A = 1, (15) 
where rj = A^. When such a transformation is found, it gives 

which can simply be turned into the Hamiltonian 



L = ^i^j - ^{A^vA),Mj, (16) 



i ^ 

We now need to diagonalize the matrix {A^vA)ij in order to obtain the vibrational eigenstates and eigenvalues. Thus, 
we need to solve the eigenvalue equation 

{B^A^vAB)^.=uf5ij, (18) 



where the matrix B must satisfy B = 1. The energies associated to the eigenstates Bij^ B2j, et Bsj are then 



= (^n, + hui. (19) 

However, in our study, we set the zero-point energy to zero, in other words we set 

Ef^ = hmuji (20) 

because the zero-point energy of vibrational modes is ill-defined. This is procedurally how we obtained our eigenstates 
(eigen-modes) and eigenvalues (eigen-energies) from our numerical solutions. But, before we discuss the results of 
our vibrational analysis oi B = 1 solitons, we must briefiy describe what exactly are the numerical solutions we are 
working with. For simplicity, we shall from hereon identify vibrational frequencies uji as vibrational energies although 
they are not exactly the same (see eqn. (fT9]) ). 

Beginning from the Lagrangian ([T|), we added the mass term ([7j), yielding the static energy 

Tr (LiV) + Tr ([Li, L,f ) 
16 V * ; ^ 32g2 VL J\ } 

.Tt(\-V -DCU^ -V^)), (21) 



8(1 -5L)) 

where i and j run over spatial components only and B is the baryonic number. The minimal energy Skyrmion for 
5 = 1 and B = 2 turns out to have spherical and axial symmetry respectively. Since we are only interested by these 
values of B, the general solution will be cast in the form of the axial ansatz 

a = i/js TTi = i/ji cos n6 7r2 = V^i sinnO tts = '02 (22) 

introduced in [l^ where il){p,z) = (V^i, '02, V^s) is a three-component unit vector that is independent of 0. The 
boundary conditions at infinity implies that '0 (0,0, 1) as ^ ^ oo. Moreover, we must impose that 0i = 
and 9^02 = 9^03 = at p = 0. 

With the axial ansatz ([22]) . we can also set the scaling length and energy in units of 2^/2 /cFt^ and Ft^/2^/2 
respectively^, giving expressions for the static energy and the baryon number that read 

Ef = - I d\Cs 



2p^ 

, , n - Wio'i (^ -I- n (^ - 

2' p2 (1-5D) 



1 ib"^ 2/3^ 1 

'a.t/; X dpxPl^ + r?^ + 7^^^ (1 - tZ-s) (1 + (1 - 2t/-3 - 4^1)) | (23) 



^ We have used 2^/2 /eFj^ and Fj^ /2^/2 as units of length and energy respectively. 
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Tl f 

B = - dzdpipi\dpip X dzip\ 



(24) 



with p = 

Also, for alfour minimizations, we fixed the constants F^^, e, and to 



= 129 MeV, e = 5.44MeV"\ and 



138 MeV. 



The values of and e were set according to ref. [8] for comparative purposes. At such values, the Skyrme model 
reproduces the mass of the nucleon and of the delta when one discounts the mass term ^ . Regarding the mass term, 
ttItt is set to its experimental value and the parameter D varies in order to probe how the mass term in ([2T]) might 
affect the Skyrmion's properties. Of course, when the mass term is switched on, the prediction for the nucleon and 
delta masses will deviate from their experimental values. This could be corrected with an appropriate choice of F^^ 
and e, however here we shall retain = 129 MeV and e = 5.44 MeV~^, since we are more interested in certain ratios 
of energies than in the actual values of energies, and comparisons with earlier works are easier. 

Physical states such as the nucleon and the deuteron require an additional contribution to the energy, the rotational 
and isorotational energy due to spin and isospin of these states. For the nucleon, this is done by fixing the quantum 
numbers of spin and isospin to / = J = ^, and assuming axial symmetry, it leads to the total isorotational and 
rotational energy of the form 



Vii 



Uii 



1 

ca7 



2f/33 



(25) 



Similarly for the deuteron one gets 



Tp2 _ 

^rot — rr ' 



(26) 



Here C/n, C/33, Vii, and Wn are moments of inertia which follows the definition in the works of Houghton and Magee 
pj| and Fortier and Marleau [l2j. Accordingly the components of these inertia tensors are 



fF, 

\2 



C/11 = 27r 



2i^l ■ 




J dzdppxPf {dpTP ■ dpTP + d,7P ■ d,7P + 2) , 

I dzdppi^lpd^i^ - zdpi^f (^1 + n^^) + z\^^ + 1{P'+ \dp'^ X d.^Pf], 
J dzdppi^[tpi {pdzip2 - zdpip2) - ip2 {pdzipi - zdptpi)] (^^ + \ {dzipaf + {dpipsf 



(27) 
(28) 
(29) 

(30) 

(31) 



+ Y (^^-^^3 + pdpipz) [dpip2dzipi - dptpid:,ip2\ + -^Y^ (2 + dptp ■ dpip + d^ip ■ d^tp) 

(32) 

Finally, the nucleon and deuteron mass are the sum of the appropriate static and rotational energies + with 
5 = 1,2 respectively. 

Having fixed the model's parameters, we then used the algorithms of simulated-annealing to accomplish the min- 
imization of ([2T]) for values of D ranging from to 0.1999. As examples of one of our minimizations, the energy 
density of the 5 = 1, 2 Skyrmions for values of the parameter = and D = 0.195 are shown in Figure [H For each 
minimization, we used a 250 by 500 point grid, representing the standard cylindrical coordinates z and p respectively, 
which provides sufficient detail for analyzing rotational and vibrational modes. Axial symmetry is therefore implied 
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FIG. 1: Energy density profiles for B = 1 and B — 2 and for values of parameter D = and D = 0.195 (top and bottom line 
respectively) on the p — z plane. Here p and z span the region from -1 to 1 in units of 2^/2/ eF^^. 



as it is known to be a symmetry of the B = 1 hedgehog and B = 2 toroidal static solutions. These configurations 
are confirmed by the results in Figure [1] along with the observations that our choice of mass term ([7j) leads to non 
shell-like energy densities. As one might have expected from larger mass terms, we also see the size of the Skyrmion 
decreases as D increases. Once a solution is reached for 5 = 1, we compute the matrices and Vij in equation 



(p!4|) numerically. We proceed with the appropriate scaling transformation, perform the diagonalization, and finally 
calculate the eigen-energies of each vibrational mode. By the form of our diagonalization (p!8|). it is clear that we would 
obtain strictly three eigenvectors depicting three types of vibration for each solution. Strikingly, but comprehensibly, 
no matter the value of the three types of vibration obtained were always of the same form. 

The first type observed may be identified to the well-known breathing mode described by the eigenvector 



(33) 



in Cartesian coordinates. The energy associated to this mode, ajtr, along with the total gives a clear picture of how 
the rigidity of the soliton changes as a function of D. Put clearly, the lower the ratio 



6r,2,3, 



(34) 



the more the soliton is expected to be deformable, or malleable. In essence, the ratio indicates whether or not 
the mass term ([7j) renders a more rigid Skyrmion, which would be desirable for the reasons mentioned earlier. The 
second type, which we call 6^2, can be understood as a vibration along x together with a simultaneous and opposite 
vibration along but not z. Its eigenvector has the form 



1 

71 



(1,-1,0). 



(35) 



The third vibrational type is characterized by a positive vibration along x and together with a simultaneous opposite 
vibration along z. We call this last type cos with eigenvector 



^(1.1.-2). 



(36) 



These eigenvectors (|33|) . (|35]) . and (|36|) . however, are idealizations of what we actually numerically obtain, even 
though our numerical approach does come close to this. For example, the three eigenvectors obtained for D = 0.8 
were (neglecting normalization) 



/ 4.064 \ 
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4.210 ^ 
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2.446 \ 


4.064 , 




-4.210 


, and 




2.446 


y 4.093 J 
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-4.869 / 
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TABLE I: B = 1 vibrational energies versus D without rotational minimization (MeV) 



D 


En 




^rot 




UJbr 


002 


0^3 





1011, 


,53 


26.67 


112, 


,12 


899.41 


290.08 


648.42 


646.63 


0.02 


1014, 


,34 


28.14 


113, 


,05 


901.29 


293.47 


651.88 


649.89 


0.04 


1017, 


,94 


29.85 


114, 


,35 


903.59 


297.78 


656.55 


654.66 


0.06 


1021, 


,31 


32.51 


114, 


,84 


906.47 


302.15 


659.20 


657.03 


0.08 


1026, 


,87 


35.38 


116, 


,64 


910.24 


308.67 


665.76 


663.10 


0.1 


1034, 


,12 


39.31 


118, 


,78 


915.34 


317.59 


674.17 


671.96 


0.12 


1044, 


,77 


44.65 


122, 


,12 


922.65 


330.16 


686.46 


684.91 


0.14 


1061, 


,48 


52.51 


127, 


,41 


934.07 


349.54 


706.13 


704.38 


0.16 


1091, 


,90 


65.36 


137, 


,33 


954.57 


384.15 


742.28 


740.63 


0.18 


1164, 


,35 


93.18 


159, 


,96 


1004.39 


463.69 


824.65 


823.08 


0.19 


1268, 


,19 


128.15 


191, 


,03 


1077.16 


573.29 


937.29 


935.27 


0.195 


1408, 


,88 


170.36 


231, 


.13 


1177.75 


716.27 


1082.78 


1078.47 


0.1975 


1591, 


,37 


220.0 


280, 


.32 


1311.05 


897.55 


1265.10 


1261.77 


0.199 


1907, 


,78 


298.78 


361, 


.49 


1546.30 


1194.68 


1567.33 


1546.77 


0.1999 


3193, 


,16 


578.49 


657, 


.96 


2535.21 


2311.62 


2726.74 


2660.59 



Throughout ah our results, an accuracy of this type was the norm. 

Table [J gives the energies of each vibrational mode, together with the static and rotational energies of the B = 1 
soliton. Here the solution was found by minimizing the static energy Eg. We first note that the energies of the 002 
and ujs modes are quite similar. Remembering that the only difference between these is a vibration along we see 
that this vibration along z diminishes the total energy of the vibration ever so slightly. Furthermore, the energy of 
the breathing mode rises much faster as a function of D than the other two vibrational modes. Lastly, we notice that 
the energy of the mass term, E^ ^ increases significantly with D. Yet, even with D = 0.1999, the mass term represents 
less than 23% of the total static energy, which is somewhat surprising if we consider the factor (1 — 5D)~^ in front of 
the mass term. Note that for D = 0.1999, we still observe a typical hedgehog configuration similar to that of Figure 
[1] although we are at the frontier {D = 0.2) of a breakdown in our numerical procedure. Indeed, the relatively large 
deviation of UJ2 with respect to ujs for D = 0.1999 could be an early sign of such a breakdown and caution is advised 
when interpreting this last set of data. 

A second set of computations were performed in parallel to the ones given above by including (iso-) rotational energy 
in the energy minimization process and therefore allowing for axial deformation of the nucleon. Such effects which 
have been shown to be significant [13, llll, [lH • Thus, adding the (iso-)rotational energy ([25|) to the static energy ([2T]) 
and minimizing the nucleon mass instead of the static energy 

En = El + E];^^ (37) 

will lead to an axially deformed solution and a new set of data shown in Table [TTl We see again that the data for 
D = 0.1999 exhibit a rather peculiar behaviour especially regarding the value uoi^r = 71.98. Although there is no 
obvious signal of a numerical breakdown in the energy and baryon number density configuration, in our opinion the 
data cannot be trusted. It should be noted that the numerical breakdown is much more drastic for higher values of D 
where it is generally characterized by scattered densities over the entire grid. The D = 0.1999 data set is nonetheless 
listed to illustrate the signature of a possible numerical breakdown. 

Notice how little the rotational energy changes with D and how similarly the vibrational energies behave as a 
function of D with those obtained without rotational term in the minimization. However, there are significant 
differences between the two sets of data which are best exhibited in a plot of the relevant ratios for both types of 
minimization. This is shown in Figure [2l 
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TABLE 11: B — 1 vibrational energies with rotational minimization (MeV) 

D En El 

990.23 51.52 73.82 916.41 269.61 541.87 540.31 

0.02 992.81 53.74 74.19 918.61 274.01 545.53 543.54 

0.04 995.98 56.42 74.65 921.33 278.81 549.15 547.23 

0.06 999.83 59.94 74.87 924.96 287.45 555.88 553.40 

0.08 1004.83 64.24 75.24 929.60 297.15 563.21 560.02 

0.1 1011.55 69.77 75.73 935.83 309.24 572.19 568.21 

0.12 1020.92 77.35 76.04 944.88 329.01 586.01 579.82 

0.14 1035.27 88.27 76.28 958.99 356.36 603.49 595.97 

0.16 1060.49 105.10 77.66 982.82 395.82 629.89 616.99 

0.18 1118.28 139.27 77.67 1040.61 493.25 680.96 664.99 

0.19 1199.51 178.84 79.68 1119.84 603.37 736.22 718.97 

0.195 1309.04 224.52 83.31 1225.72 729.46 803.61 783.27 

0.1975 1452.83 277.16 89.70 1363.13 808.13 904.12 895.73 

0.199 1707.07 360.53 103.82 1603.25 1003.54 1135.50 1043.64 

0.1999 2309.34 578.49 151.63 2157.71 71.98 2102.43 1620.67 
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a) b) 

FIG. 2: Ratios of vibrational with respect to the static energies Ri = ^ as a function of the mass parameter D where in a) 
i labels the vibrational mode i = br and in b) i = 2. Triangles (squares) correspond to the solutions without (with) rotational 
energy in the energy minimization. The vertical solid line in a) is the experimental value of ratio — where Enoeper is the 
energy of the Roeper resonance, sometimes identified with the breathing mode of the proton, but shown here for comparison 
purposes only. 



9 



The differences between the ratios of the minimizations with and without rotational energy are now clear. For 
the breathing mode, the curves overlap quite well. However, for the other two types of vibration, UJ2 and c^s, the added 
rotational energy to the Lagrangian lowers significantly their energies. Perhaps because E^^* renders a slightly less 
oblate Skyrmion, the vibrational energies along p diminish accordingly. The general tendencies of all curves however 
indicate that the mass term does add a certain "rigidity" to the field configurations as D increases. What remains to 
be seen is if there is a correlation between this defined rigidity and the binding energies of B = 2 solutions. This is 
done in the next section. 

But before analyzing binding energies, we must mention here the complementary results of Lin and Piette [I^. 
The authors consider the vibrational modes of the Skyrme model with the mass term falling into the class of eqs. ([5]) 
and namely of the form Tr{U^ — 1) which differs from ([7|). They follow the time dependent approach introduced 
by Barnes et al. [17] to identify the vibrational modes which assumes local instead of global scale transformations 
that we chose to perform our calculation. Also they use of the rational map ansatz which seems to prevent a precise 
determination of the vibrational energies. For anyone of these reasons, their results can not be compared directly 
with ours. Nonetheless they observe a clear increase in the vibrational energies with respect to the mass term which 
is in qualitative accord with our results. 



IV. BINDING ENERGIES 



The results for the static energies of the B = 2 solitons are presented in Table Hill together with their corresponding 
binding energies computed from the total static energies of the B = 1 solitons of Tables [H Also shown are the results 



TABLE III: B — 2 energies without rotational minimization (MeV) 



D 


Ei 


2Ef 


-El 


Ed 


2En - 


-Ed 





1720.55 


78 


27 


1813.51 


209 


55 


0.02 


1724.91 


77 


66 


1818.19 


210 


49 


0.04 


1730.18 


77 


01 


1824.65 


211 


23 


0.06 


1736.7 8 


76 


15 


1832.75 


209 


86 


0.08 
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for the physical states of deuteron Ed versus that of the nucleon En- All results in Table HTTl were obtained by 
minimizing static energies alone, i.e. minimization without rotational energy, to avoid non-axial solutions. We recall 
here that we have fixed the F^^, e and m^r parameters of the Skyrme model leading to results for Ed and En that 
are much higher than their experimental values. While an appropriate fit of the parameters could fix this problem it 
is not necessary in the context of this work since we are more interested in the relative weight of the vibrational and 
binding energies than in their actual values. With the data provided by Table HIH we define i?, the ratio of the energy 
of the B = 2 soliton with respect to that of two isolated B = 1 solitons. R is plotted in Fig. [3] for three sets of points. 
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corresponding to (iso-)spinless solitons and 
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Ed 
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(39) 
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for comparison of the deuteron versus two nucleons. These quantities indicates how the relative binding energies 
change with increasing D. The dashed hue represents the experiment ahy measured mass of the deuteron over twice 
the mass of an individual nucleon {R = 0.9998). 

Clearly, as D approaches its critical value of 0.2, we observe that the binding energies for / = J = solitons 
(triangles in Fig. [3|) diminish considerably without ever crossing the dashed line or the value R = 1 which corresponds 
to the limit of instability of the B = 2 solution. The sharp increase in R suggest that it may be possible to adjust 
the value of D such that the relative importance of the binding energy would be arbitrarily small — of course, the 
solution would still be of toroidal form which is presumably not that of the deuteron. Some cautionary comments are 
in order here. First the last point for this set of data shows a sharp increase of the binding energy Sit D = 0.1999. 
It is at the boundary of the region of breakdown of our numerical technique and attempts to obtain stable solutions 
as we push D closer to 0.2 were not successful. Also the limit D 0.2 is ill-defined and approaching this limit is 
physically questionable as the contributions of the mass term which breaks the SU (2) symmetry gets relatively large. 

Surprisingly, the situation is more obscure for the deuteron-nucleon data. A first set of data, based on static energy 
minimization (squares in Fig. [3]), shows that the relative importance of the binding energy is quite insensitive to D 
up to 0.17 but increases (i.e. R decreases) dramatically as D approaches the critical value of 0.2 contrarily to what is 
observed for the spinless Skyrmion. This is entirely due to the distinct behavior of the B = 2 and B = 1 rotational 
energies. Of course rotational deformations may affect this behaviour. To provide some measure of that effect we 
present a third set of points on Fig. [3] (circles) where R is computed with the B = 2 static and B = 1 rotationally 
deformed solutions. Note that an exact full 3D B = 2 rotationally deformed computation would lower R^ and so this 
third set of data represent the absolute maximum of R for the deuteron-nucleon system. Again the mass term ([7j) 
that we analyzed always lead to bound states {R <1) which is an interesting result by itself. However the observation 
of this opposite behavior with respect to the parameter D for these last two sets of data emphasizes the importance 
of the rotational deformations. It may well be that a full 3D computation oi B = 2 deformed solution leads to a 
behaviour closer to that of the second set of data (squares) but our results does not allow to infer on the exact nature 
of the bound state for the deuteron. Perhaps a full 3D computation allowing for non-axial deformation of the deuteron 
would conclude otherwise. It is also clear that the bound states of spinless Sakyrmions and deuteron-nucleon system 
may show completely different behaviour. 



0.05 



0.1 0.15 0.2 



FIG. 3: Ratios 2^ (triangles) and (squares) as a function of the mass parameter D. The dashed line at = 1 corresponds 
to the limit of instability of the bound state whereas the experimental value for = 0.99979, i.e. very close to 1. 

Additionally, these results put into question the conclusions obtained by Piette and Zakrzewski Q. Their results 
show that for a D larger than approximately 0.12, the toroidal configuration is an unstable one because its total 
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energy is larger than twice that of a single nucleon. More precisely, their results indicate that the ratio crosses 1 
ai D ^ 0.12 and continues to increase afterwards, in contradiction with our results (Fig. [3]). This discrepancy suggests 
that the use of the rational map ansatz does not adequately minimize the energies of the B = 2 configurations, and 
that we cannot have too much confidence in its ability to truly find the minimal solutions of any B > 2. Hence, the 
more exact simulated-annealing approach for finding solutions may be the best way to decipher which model types 
are more suitable than others. It may also be the only method for finding out why the B = 2 Skyrmion has a toroidal 
shape in the first place. There is also certainly a link to make between the increasing rigidity (characterized by ^) 
and the changing binding energy of the solitons as a function of D which also depends on the nature of the solitons, 
(iso-)spinless or deuteron-nucleons system. In our view, both characteristics are needed for a better model, and hence 
adds to the appeal of further dynamical or mass terms such as ([7]). 



V. CONCLUSION 



From the behaviors of the binding energies (Fig. [3]) and the breathing mode energies (Fig. [2^), we can assert that 
the mass term ([7]) with a non-zero value of D does not quantitatively succeed in providing a consistent model of 
nuclear matter in the sense that the experimental value for binding energy of the deuteron was not attainable for the 
values of D considered here. It also emphasize the need for rotationally deformed B = 2 computations to determine 
the exact behaviour as R seems to be very sensitive to rotational deformation near D = 0.2. Meanwhile it may be 
premature to conjecture the exact deuteron-nucleon behaviour on the basis of the results for spinless solitons. Our 
results have also shown that in contrast to the rational map solutions of Piette and Zakrzewski, the behavior of the 
ratio 2^ never crosses 1 and seems only to approach 1 in the limit of D — > 0.2. Hence, this indicates not only that 
the 5 = 2 is a bound state for all values of D considered in this work, but it also demonstrates the clear limitations of 
the rational map ansatz in providing exact quantitative insight into how the model behaves and works. In conjunction 
with the behavior of the binding energies, the ratio of ^ shows that as D increases the field configurations are more 
resistant to being vibrationally excited, which is indicative of more rigid solitons and perhaps a relation between 
rigidity and larger binding energies for the deuteron. 

All calculations relied on fixed values for the parameters of the Skyrme model (except for D) mostly for comparison 
purposes with previous work. These parameters are usually fitted to reproduce the experimental values of the mass of 
the nucleon and delta or other physical quantities. Indeed, the values of the energies we have obtained were sometimes 
far from those of experimental data. Therefore, some of our conclusions may no longer hold for a more physical choice 
of Skyrme parameters despite the fact that they were based on the relative importance of each quantity. On the other 
hand, three dimensional simulated- annealing programs will be needed to explore baryonic numbers beyond 2, and to 
confirm the validity of axially symmetric solutions for B = 2. Regarding the prospect of obtaining 2-nucleon shape 
solitons for 5 = 2, it would seem that the mass term in ([7j) is not sufficient and perhaps we need to rethink what type 
of lagrangian would permit such a minimal configuration. In any event, it is also possible to add as many parameters 
in the mass term ([7]) with higher powers of U as we wish, and this alone may lead us towards a sounder effective 
theory of QCD. 

This work was supported by the National Science and Engineering Research Council of Canada. 
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